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TECHNICAL MEMORANDUM 14-32 


EVAPORATION, HEAT TRANSFER, AND VELOCITY DISTRIBUTION 
IN TWO-DIMENSIONAL AND ROTATIONALLY SYMMETRICAL 
LAMINAR BOUNDARY-LAYER FLOW* 

By Nils Fros sling 

INTRODUCTION 


Aside from the simple case of the plane, no quantitative calcula- 
tions of the evaporation of a "body in a moving medium exist so far. The 
heat transfer, which under certain circumstances (see p. 4) follows the 
same laws, has been treated theoretically for the cylinder by Kroujillne 
(ref. 1) and Squire (ref. 2). For boundary-layer flow, Kroujiline used 
for the temperature field a power-series method of the same type which 
has been Introduced for velocities by Pohlhausen (ref. 3 ) • Because of 
this stipulation of the profile form, the result must be approximate, 
and the eventual agreement with the correct value is rather accidental. 
Squire gave an exact treatment of the transfer in the immediate proximity 
of the stagnation point. However, it is of great interest to have a 
calculation method for heat and mass transfer in the entire boundary 
layer, the error of which depends only on the work expenditure of the 
numerical calculation and, therefore, not on possible approximative 
formulations. Even though the calculation is time consuming, one has 
the advantage of being able to check approximate and more rapid methods 
with respect to this solution. Under the supposition that the constants 
of the problem (shape of body, pressure distribution, etc.) may possibly 
be eliminated from the equations to be solved. It is also possible in 
several cases to -use the complete exact solution directly. The author 
of this report perfected, for this reason, the exact solutions for the 
temperature and concentration fields. Two dimensional and rotationally 
symmetrical steady boundary-layer flows were treated. The latter case 
is the more complicated one because of the form of the continuity 
equation. 


* "Verduns tung , Warmeubergang und Geschwlndigkeitsverteilung bei 
zwe idimens ionaler und rotations gymmetrischer laminar er Qrenzschicht- 
strdmung. " Lunds Universitets Arsskrift, N.F. Avd. 2, Bd. 36, Nr. 4- 
Kungl. Fysiografiska Sallskapets Handlingar, N.F. Bd. 51, Nr. 4, 194-0. 
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For the calculation of the transfers, however, the velocity fields 
must be known. For the two-dimensional case one has, aside from approx- 
imate methods of solution by Pohlhausen (ref. 3), K^rm^n and Mill ikan 
(ref. 4), and others, the method of Blasius (ref. 3) and Hiemenz (ref. 6) 
which was Improved by Howarth (ref. J). By means of this method one may 
solve the equations, without any arbitrary assumptions regarding velocity 
profile and the like, by power-series development from the stagnation 
point up to an arbitrary point on the meridian curve. The work expend- 
iture depends on the required accuracy and on the position of that point. 
Since the development becomes very rapidly more cumbersome with the 
distance from the stagnation point, it is appropriate to use, from a 
certain point onward, a continuation method, for instance, according to 
Erandtl (ref. 8) and Gortler (ref. 9) • Howarth (ref. 7) transformed 
the functions of Blasius and Hiemenz into functions of such a type that 
the constants disappear from the equations so that the numerical solu- 
tions for them may be applied to any two-dimensional flow. He treated 
the symmetrical as well as the unsymmetrical case, and Indicates the 
solution of one of these functions. For the present investigation, 
Howarth' s functions are used in the two-dimensional case. Since the 
accuracy of Howarth 1 s numerical tables is not sufficient for the cal- 
culation of the transfer, a new numerical calculation was made of certain 
functions. For the three-dimensional rotationally symmetrical case, in 
contrast, there exists, so far, no calculation of the functions of the 
series development. Because of the modified continuity equation, other 
systems of equations must be used, and these systems are established 
here. The necessary functions are numerically calculated. Although the 
form of the meridian curve takes effect through the continuity equation, 
one can proceed in the distribution of the functions in such a manner 
that the constants of the meridian equation disappear, and the solutions 
therefore are valid not only for arbitrary pressure distribution but also 
for arbitrary shape of the body of revolution. The continuation method, 
beginning with the limit of validity of the broken-off power series, has 
been perfected also for this case. 

An investigation is carried out regarding the validity of the law, 
stated, for instance, by Ulsamer (ref. 10) that the Nusselt number is 
proportional to the cube root of the Rrandtl number. 

A few approximation methods for the calculation of the transfer 
layer are discussed. 

Only a brief survey is presented here since a more detailed report 
is to be given in a later paper. 
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THE FUNDAMENTAL EQUATIONS 


With dimensioned quantities, the boundary-layer equations for flow, 
concentration, and temperature read 


3u 

U iiii + v --H = ug' + V 
bx by by* 


— + — = 0 two dimensional 
3x by 


B(ur) S(vr) 

+ = 0 rotationally symmetrical 


dx 


Sy 


2 

3c , be . c 

U + T = A — S' 

3x by by 2 


u + v = a 
bx by by 2 


Concerning the derivation of the two last equations for the 
rotationally symmetrical case see p. 15-16. The boundary conditions are 


y - Oj u = v = Oj t = t 0 J 0 = 0 *} 
y = a>} u = U} t = 0} c = 0. 


In these equations the customary designations of the various 
quantities are used - , x = distance along the body surface from the 
stagnation point to the base point of the normal to the body surface, 
y = length of that normal, r = distance from base point to axis of 
rotation, u, v = velocity components in the direction of x and y, 
respectively. U = the velocity component parallel to the body surface 
immediately outside of the boundary layer (calculated from the experi- 
mentally determined pressure distribution), t = excess of the temper- 
ature of the surface over the temperature of the undisturbed fluid, 
c = corresponding concentration quantity, v = kinematic viscosity, 

a = temperature diffusivity 



A = diffusion coefficient. If 
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Uq is the undisturbed velocity and D a characteristic length of the 

body (for instance, its diameter), one can transform the equation Into 
dimensionless form by dividing the velocities, lengths, temperatures, 
and concentrations by the quantities Uo, 33, to, and c m . These equa- 
tions contain as constants, among others, the Reynolds number Re = 

This number one may eliminate by modifying the scale of the boundary 
layer in transverse direction, by multiplying the values of y and ^ 

by x/ie . The equations used below with the designation "dimensionless 
equations without Reynolds numbers" are changed in their appearance, 
compared to those mentioned above, only in that v disappears, and A 

and a are replaced by A and &, in the boundary conditions to and 

A 

c m are replaced by 1. The two quantities — and — which are often 

independent of pressure and temperature, as in the caBe of ideal gases, 
are dependent on the media used. These quantities are called Stanton's 
numbers. Frequently their inverse values are used, designated as Rrandtl 
numbers . In an earlier report of the author on the evaporation of drops 
(ref. 11) the designation a is used for the Stanton number. Since at 
present this letter is used mostly for the Rrandtl numbers, this defini- 
tion is employed in the present report to prevent misunderstandings. 

Thus a here signifies: a = jr or, respectively, = ^ = — 

A a \ 

The equations for temperature and concentration are therefore 
identical when t and c are interchanged. For the temperature- 
boundary layer it is assumed, however, that the dissipation and the 
heat generated. by change In pressure may be neglected. This assumption, 
is satisfied for not-too-large velocities (ref. 12). The equations also 
presuppose that the velocities be small oompared to sonic velocity in 
order to make the compressibility negligible. A further limitation of 
the equations is given by the fact that the differences in concentration 
and temperature must not be so large that the constant characteristics 
of the media vary from point to point. Because of the identical form 
of the two equations for temperature and concentration, which is thus 
satisfied under these presuppositions, both may be treated simultaneously. 
In the following equations one may, therefore, immediately interchange 
the quantities c and t. | 

In the search for a solution which satisfies the accuracy require- 
ments discussed in the introduction, the method of power-series develop- 
ment In x was used. Breaking off the power series after a certain 
number of terms one was able to use this solution from the stagnation 
point up to a point the position of which was dependent on the accuracy 
requirements. Starting from this point one could then use for the layers 
of different types step-by-step continuation methods. 
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For the sake of brevity, ve shall use below the common name 
"transfer" boundary layer for the temperature and for the concentration- 
boundary layer. 


POWER-SERIES DEVELOPMENTS TS x 


A. TWO DIMENSIOHAL CASE 
a. Flow Boundary Layer 
1* Symmetrical case 


As was mentioned In the introduction, this case has been treated 
by Blasius, Hiemenz, and Howarth. This report uses for the most part 
the same designations as Howarth. The only difference is that Howarth* s 
quantities F v are here replaced by the quantities because the 

capital letters are more suitable for the functions of the transfer 
boundary layer. 

In order to replace the two unknown quantities u and v by a 
single one (t ) , the following conditions satisfying the continuity 
equation are set up as usual: 


U Sy V 3x 

The first flow equation then becomes (with dimensionless quantities 
without Reynolds numbers) 

St cfi = uu, + 

by bxby dx $y 2 dy^ 

Blasius and Hiemenz solved this equation by means of the formula 


f = f-jX + + . • . 

for the symmetrical case where the velocity distribution outside of the 
boundary layer follows the formula U = ujx + u^x^ + u^x^ + . . . . t v 
are functions only of y. 
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R y comparison of the various powers of x the equations for i|r v 
were obtained. These equations were freed of the constants Uy by- 
introduction of the functions f v , g v , etc., by means of the following 
statements : 


*5 = 


__5 

nFi 


Tl = 7^ 
, g 5 + ^ ^ 


10u, 


^9 = 


>/* 


l g 9 


U5U7 

U1U9 




*3=^ f 3 


+7 


= 8u T t 




UzUc yx-z/ 

te7 + ^ h 7 + ^ k T J 


b.9 + 


u 5 2 v u ?S 
^9 + 


U7 


u i u 9 


u 1 -u 9 


J9 


u-^uo 


*9 


2. Uasymmetrical case 

For an unsymmetrical two-dimensional body for which the velocity 
distribution follows the formula U = u^x + Ugx^ + u^x3 + . . . , the 
formulation 


t 


= i|r-j_X + 'if 2^ 



+ . . • 


was used. Here also the were freed of the constants Uy, this time 

by the expressions 


n 



tl = f lv/^i 



etc. 


b. Transfer Boundary Layer 
1. Symmetrical case 

The author of this report attempted in the dimensionless equations 
without Reynolds numbers, aside from the formulations for ijr mentioned 
above, a development for c in the following manner (c v are functions 

of y only) : 


-£ 

v=0 


C v xV = Cq + C3X + C 2 X 


•2 + o 3 x’ 


CljX^ 


c 
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Boundary conditions: 

y = 0 cq = 1 = Cg = • • . = 0 

|y=co cq=c^ = C 2 = ...=0 


Ey substitution one obtains for the c ^ and c 2 q+± 


n n 

” c "2n = 22^ ' 2n+l-2k c 2k - 22 (2n + i - 2k)t2n +1 >2k c *2k 


- c" 

a 2n+l 


n 

>_ 

k=0 


- C 22 C 2 ^ + 1 ^ , 2n+l-2k c 2k+l 


(2n + 1 - 2k)t 


2n+l-2k C ' 2k+l I 


One can easily show that the equations for c 2n+l axe suc ^ that 

they become identically zero. In the groups of equations mentioned 
above -which constitute the recursion equations, there occur exclusively 
functions with even or odd subscripts i From this one can see that c 
is an even function of x which follows, besides, from the nature of 
the problem. In order to be free of the constants TJy, new functions 

are introduced by the following statements: 


*1 = y \Ju±} + v as 


above: 


c 0 


It-UxFo 

F 0 ’ °2 - 


_5.| 

u i 


% 


U* 


8 


u l«5 


Ut \ 6 u-jju^ 


=6 


u*3 


U l 2u 7 


*6 


10u^ 


l 


u^ 


Us 


u ?% 


«5 


C 8 = ~\ G 8 + 35 x 5 =8 + ^*8 + ^ J 8 + ^ y 


Boundary conditions: 

tj = Oj Fq = 1; the remaining functions = 0 


n = °°; 


all functions = 0 
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h" 

- -f x f 2 

' + 2fl'F 2 

" » 3 f 0 

t 





k" 

= - f iV 

+ ^ 1 % 

- ^Fq 1 

1 





k" 

- -w 

+ bt^ 

- 5h 5 V 

** 

(2fj’F £ 

- 3 f,v)_ 



k" 

- -f 1G6 ' 

+ 6 f 1 'G 6 

“ Tg 7 F 0 ’ 






k" 

= “ f l H 6 ' 

+ 6f l' E 6 

-TVo’ 

+ 3 ( 

'W 

+ 2g 5 'F 2 

' 5 W 

- Js^) 5 

k" 

= - f l K 6’ 

+ 6f *K 

1 o 

“ 7 Vo' 

+ 3 ( 

4f 3 % 

+ 2h 5 'F 2 

VW 

' Va)j 

lr« II 
0^8 

= -w 

+ 8f 1 , 0g ■ 

‘ 9 BgV 








+ • 

■ sw 

+ f^i ' a 6 

+ 2g ? 'F 2 

- 5 W 

- 76 f 2 ') 

k>"- 

= -w 

+ - 

■ Vo* 

+ # 


- 



2£ 

■i 

‘ “ f 1 J 8 ' 

+ - 

• 9 W 


( 6f 3 ,H 6 

+ 2h 7 'F 2 

' 5 W 

- 7 b/ 2 ') 


f( h 5 ,Q 4 + g 5 '^) 

- 18(^5^ • + 

• «*) 




II 

& 

H| b 

' -f lQ8 - ■ 

- ^'*8 - 

9W 

*¥{ 

h ,K 6 

+ 2^'Fg 

- 5 W 

- 7V 2 ') 


fi kh 5'\ - »sV) 


O^ese functions are simpler than those for the flow boundary layer 
insofar as the equations are only of the second order. gsJef iSt eaua 
tion may be solved by quadratures: trrst e<lua " 
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*0 



The re m aining equations do not have explicit solutions and must therefore 
be integrated by other methods, for instance, according to Runge and 
Kutta. (See p. 19.) 


2. ^asymmetrical case 

With c = Cq + c^x + CgX^ + Cjx3 + . . . and the boundary conditions 


y = Oj Cq = lj c i = c 2 = • . • * 0 

y = 00 . c 0 = c-^ = c 2 = . • . = 0 


one obtains for c v 



v 

y~ kfr ' c 

fco v+1 ~ k k 


t 


(v + 1 - h)^. 


v+l-k c k 


Here none of the functions c v disappears, and c is therefore, 

as had been expected, an even function. The equations here are not 
divided into two independent groups but the functions follow successively 
one from the other. 


Distribution of the c v : 


n = y^J +1 = f iv/^ etc -J c o = fqj c i = 




ug 


5% 


"2*3. 


02 = — — IG2 + - B g/j Cx = 1G3 + H* + 

U 1 \ u l u 3 )■ ■> u i V ^ u l% ” u^u^ 


^2 1 ; 
ui^ 7 




= 5^1 
U 1 


112% 


U 3 


V5 V 5 


% + + + » 


u 2 


u l ?u 5 
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Boundary conditions: 

il = 0; Fo = lj the remaining functions = 0: 

< 9 

T = «; all functions = 0; 

5 F 0 " - -Vo' 

5*1 ” ‘ - f i F i ' + fi' F i - 2f 2 r 0 1 

' “W + * 1 % - 3« 3 V 1 

= ' f lV + - 3 Vo' + f(V F l - 2f 2 F 1 ')| 

&5" . -fjOj' + 3fyG, - IJajFo’ & 

0 H 3 = -f 1 H 3 1 + 3f 1 'H 3 - 4h^F 0 ' + ^(2f 2 'G 2 + g 3 'F x - 2f 2 G 2 ' - 3 g 3 P 1 ^ , 

■ ’W + * 1 '^ - Vo' + ^(^ 2=2 + V F 1 - - 3hj F i) j 

K" ■ - f iV + 4 f i'\ - 5g 5 P 0 ' 

■ -W + - 5h 5 F 0 ' + |(3f 2 'Gj - 2f 2 G,' - 

h" - - f lV + - 5^ F o’ + |(2 Sj'» 2 - 3gjG 2 ') 

yH" - -V 4 ' + - 5J 5 P 0 ' + 'Hj + h lf 'F 1 - 2 

I y( S 3 ,E 2 + h 3 ’ G ^ “ 8 (* 3 * 2 ' + h 3 G 2 j 

5^" - -W + - 5^0 ■ + |(3f 2 'K 3 + . ayy . 4 Vi .j + 

case. 5516 flrSt 6qUation is ide ntical with the first one , of the symmetrical 
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B. ROTATIOHALLY SYMMETRIC CASE 
a. Plow Boundary Layer 


For flows about a blunt body of revolution whose axis lieB in the 
direction of the flow, the flow equations are, according to Boltze 
(ref. 13), 


u— + v— ■= UU’ + 
dx by 

Mur) + Mvr) = Q 

Sx by 

In the transformation to the dimensionless form without Reynolds 
number v disappears. The quantity r then must have, for the bodies 
of revolution, the meaning, distance of axis of rotation up to the base 
point of the normal instead of up to the point (x, y) . A function for 
identical solution of the continuity equation is desired. Boltze 
(ref. 13) suggested a function i|r which is defined as follows: 

u . iit T ..ih 

r 3y r bx 


In the present report another solution 4" also has been examined. 
Definition 

_ &t|r 1 d(i|fr) _ bif ip dr 1 d(\f r) 

dy r by bx r dx r &x 


The function ^ has the advantage that the equations become some- 
what simpler and that the velocity profile is obtained directly. The 
functions are derived in both cases. i}r and ijrr are flow functions. 



1. Use of the function ijr 


After substitution of the expressions for u and v into the 
first boundary-layer equation one obtains 


dA 2 dr + StJt dAlr 

hyl ax by bxby 


t & 5% = Pmj. + r 2 

by 2 by$ 
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The power-series developments used ^i|r v function of yj r v and 
u v constants ) 


ilr = + + Ifgx^ + . . . 

r = rjjx + r^x^ + r^x^ + . . . 

U = u^x + u^x^ + u^x^ + . . . 

The functions i|r v have the following boundary conditions: 

y - oj t v = V ■ °; 

y = 00} Tjfg' ■ r^j = r L u 5 + rju-L} %' = r^uj + rjUj + r^j . . 


After substitution of the power expressions into the equation for 
i|r one obtains equations for by comparison of the different coef- 

ficients. These equations may be freed of the letters r v and 
by the following formulations: 


riui 2r, u* 




r-,u 


13 


=h 4 


% * 


s 6 


r n u = ' UjU^ 


15 


r l u 5 


rx 2 Un 


3rnU 5 / rcu-, Mif- r^U7 ryu x 

+ + — k 6 + + 


• i * 


The new functions have the boundary conditions 


t) = 0; all functions and their first derivative = 0 

*1 = “i *2* = Oj.' - 1%* = 66* = V = J6* = 3 j kg' = qg' = 0}. 
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One obtains the following equations: 

f 2 ,u = -fat 2 " + |(f 2 ’ 2 - 1 ) 

V" = - f 2 %" + 2*2 V “ ** 2 % - 1 

V" - -W + “s'V - * 2 \ - i^ 2 2 - + x ) 

- 86 "’ - "W * * 2 %' - 3 V 8 S - 1 

V" - 'W + 3f 2 'V - 3f 2 "h6 - |(»2' 2 - 2fs<2” + X ) 

V" = -f 2 k g " + 3f a 'V - 5V*g + |(3^' 2 - W) - I 

J 6 "' - -W + »2'J6’ - 3 f a"J6 + HH' - f^V 4 \%”) + 

jfak" - 4f 2'%' + &2\ - x ) 

qg m - -f 2<36 ” + 3f 2 'qg' - 3f 2 ”qg + f(3^- 2 - SV’) + f^yv* - ‘MfaH 

“ 2 % + K 5f 2' 2 ■ “a'a") 

2. Use of the f met ion t 

One obtains 

dfr d 2 ^ dfr 5 s ? 1 dr - d 2 * _ , d^jF 

dy dxdy dx £y 2 r dx ^ ^y2 1111 + ^3 

Power series developments 

if - ijr-jX + ij 'r-^a? + + . . . 

3 5 

r = r^x + r^x-^ + r^x^ + . . . 

U = u-jX + u^x^ + u^x 5 


+ . . 



MCA 3M 1452 


14 


Boundary conditions 

y = 0; t v = t v ' = 0; 

y = oo; t v ' = Uy 

The functions ip v are here divided up as follows: 




U 1 




f l? ?3 = 





^ = 


5ub 


f 5 “ 


5u 5 


rcu- 


5 U 1 


s/^i *1*5 


h* + 


u r 

u l u 5 


r 5 u 5 

r l u 5 


1 + 2^ 
5 r;L 2u 5 


Boundary conditions 


tj = 0} all functions and their first derivative = 0; 
n = “i *i' - 1 ; S3' = 85' - y - V “ V ■ ^ 5 * = V = 0j 

f l"’ - “W + |( f l' 2 " l) 

83 m = - txS 3 " + 2 f 1 l S 3 l - 2f 1 "g 3 - 1 

V” = "W + & 1 'V “ - 2 

85"’ - -*185" + 5 f 1 , g 5 ’ - ?fi"g 5 - 1 

V" = - f i h 5 " + 3f l'V - * 1 *? “ I f l f l" 

k 5 '» = -f^" + 3f 1 , k 5 l - 3f 1 "k 5 + 2g 3 * 2 - | g 5 g 5 " - | 

< 35 = - f i« 35 " + 3 * 1 ' ^5' ~ 3 *i"J 5 + ^83 *h 3 ' - | (63*13" + h.363’^ - 

f (fl '*3 + f l g 3 ") 

05 ,M = + 3 ^' <35 r - + 2h 5 * 2 - | h^" + 1 f^” - 

f ( f lV + t ^) 
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la the first method, the functions have even subscripts, in the 
second odd ones. A simple relation exists between the two groups of 
functions which one may easily obtain by equating the two expressions 
defining u_ and, respectively, the two expressions defining v (which 
gives i|r = -i|nr ) . 


f 2 = f 1? 

fl 

% = e y ^ = h 3 + Y’’ 

£ 

&6 = h 6 = 115 + *6 = k 5; i 6 - h + | 63* <36 - ^5 + | *3 


b. Transfer Boundary Layer 


The general equation of the temperature and concentration fields 
for rotationally symmetrical flow has not been set up before. For the 
special case where the body is a sphere, the author (ref. 11) has shown 
that for boundary-layer flow the equation is identical with the one for 
the two-dimensional case, at least for points which do not lie directly 
at the stagnation point. In the present report, it is shown that the 
same boundary-layer equation is valid also for arbitrary blunt bodies 
of revolution, and that this applies to points directly at the stagnation 
point as well. The introduction of mass or heat into a volume element 
by diffusion and convection is expressed by the following equation (which 
is valid for rotationally symmetrical flow without neglect of the boundary 
layer when x and r are counted up to the element instead of to the 
base point) : 


The derivation becomes the simplest if one chooses as the volume element 
an element bounded by two meridian planes, two surfaces x = constant, 
and two surfaces y = constant. In order to arrive at the boundary-layer 
equation, one groups the derivatives 


2. d(ur) + c S(vr) + 
r Sx r Sy 


u — + v 


Sc 

by 



Sc 1 Sr t 
Sy r Sy 


+ 1 <k. £ 2 . 

r <3x Sx 
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Supposing that a thin boundary layer exists, the terms 2 and 3 in 
the parenthesis disappear. The first two terms of the left side disappear 
because of the appearance of the continuity equation. The last term of 

the parenthesis becomes infinite at the stagnation point if ^2- is not 

ox 

here zero. In order to avoid discontinuities at the stagnation point, 

one must therefore equate there ^2. = 0. Then the last term becomes 

ox 

everywhere negligible, and one obtains the equation 

u + v ^2. = A 
ox dy 



which thus is identical with the one in the case of two-dimensional flow. 


1. Use of the function t 


Because of symmetry and of the requirement 

c the expression c = cq + cgx 2 + c^x^ + . . . 
boundary conditions for Cyj 


^2. = 0 one uses for 
Ox 

with the following 


y = 0; c 0 = 1; e 2 = c^ = . . . = 0; 
y = ooj c 0 = c 2 = cjj. = . . . = 0 


For the c v one obtains equations which contain r v and Uy. In 

order to eliminate these constants, one may make the following 
substitutions : 


r 3 u l \ 
Go + Ho| 

r l u 3 ^ 


1 ^3 / 

= yJ 2UJ.; i v = as before} Cq = Fq} c 2 = — - G 
v u x 1 

3U5 3u= / rcUn u*^ riu* r^^u-i 

04 = — Fit = — G4 + K1 + + -^4 qU 

U 1 U 1 \ r l"5 u l^ r l^ r l‘S 


Boundary condition 


H = 0} Fo = 1} remaining functions = 0} 
T] = «} all functions = 0 
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For the new functions the following equations of the second order 
are obtained: 

i F 0 " - -t^ 0 ’ 

1 Oa” = -f 2 0 2 ' + f 2 'G 2 - 2g k F 0 ' 

i H 2” - ‘W * f 2 ' H 2 - Vo' + 2 Vo 

i V - - f 2 ^' + 2f 2 '^ - 3B6V 

1 V - -f 2 V + 2f 2 '^ - Sbgl,,' + \ 1 ^ 0 - 

i v . -fgiv + 2 f 2 'K i| . - 3*gV + y V G 2 - 1 %V 

i V = -f 2 V f 2f 2 'J^ - 3j fi F 0 ' + |( % 'S 2 + V G 2 ) - 

t(W + W) ' f (- f 2<V + f 2'°2 - Vo') 

i V = -* 2 %' + 2t 2 '% - 3^' + | h^'Hg - § h,^' - 
|(- f 2 S 2 ' +f 2 ,H 2 * Vo' + 5Vo') 

2. Use of the function 

With the same power development for c and the same def init ions 
of the functions Fq, etc., one obtains 

1 F 0 " = -f^o' 

i 02" - -fjCfe' + f]_'G 2 - 2gjF 0 ' 

i %" - -qV + q-ife - 2bjF 0 ' - | f 1?0 ' 



1 

a 


V 






-* 1 %' + - 3g 5 F 0 ' 

" f lV + 2^ -3h 5 V -fVo* 

" f lV + + | «3 ,G 2 - f g 3 G 2 ( 

- f lV + *1% ~ 3 b F o' + |( g 3 ,H 2 + h 3 ,G 2)- f( f lV + 

g 3 F 0 ') " j( g 3 H 2 ' + VV) 

~ f lV + * 1 % ' + | V H 2 “ f W “ f ( f lV + h 3 F o) + 


One can show easily by application of the relationships between the 
functions with even and the functions with odd subscripts that the systems 
of equations for the cases 1 end 2 are identical. 

The first equations of the two systems are identical with the first 
equation for two-dimensional flow and are, therefore, also solved by 
quadratures . 


C. FINAL EXPRESSIONS FOR THE TRANSFER 


Nu = 


The transfers are made dimensionless by the Nusselt number 
D d 2 !!! D B 2 Q 


Ac m BsBt 


or, respectively. 


Ato BsBt’ 


It is easily shown that 


- —== ■ = /- — j ^3-, where c and y are "dimensionless and without 
vRe \ BtyO By 

Reynolds number." The heat transferred by radiation is, of course, not 
contained in this expression. For two-dimensional symmetrical bodies one 
obtains 
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Corresponding expressions are obtained in the other caseB. As one 
can see from the equations, one may easily calculate the Nusselt number 
for arbitrary pressure distributions and body shapes which agree with 
the formulations, if one has made a numerical calculation of the func- 
tions . Unfortunately, the quantity a is left over and one must there- 
fore make different solutions for different media. As is shown in a 
section below, however, one can free the equations of a, too, if a 
is large. 


D. NUMERICAL CALCULATIONS 


For the two-dimensional symmetrical case and for rotationally 
symmetrical bodies the author has numerically calculated various func- 
tions, corresponding to the three first terms of the power-series devel- 
opments in x. The method of Runge and Kutta (ref. lk) was used for 
this purpose. This method is rather time consuming but one has good 
possibilities of determining the errors. The first function fq of 
the two-dimensional case has been calculated by Hiemenz and Howarth with 
an accuracy sufficient for this investigation; Howarth' s values are 
directly used here. For f q in the rotationally symmetrical case there 

exists a table by Hartree (ref. 15) which was set up by using a mechanical 
differential analyzer. The accuracy is here not sufficient and the first 
two derivatives also are required; for this reason the function is cal- 
culated here anew. Since the equation for fq is not linear and one 

therefore cannot find the solution by combining two particular solutions, 
it was valuable to have approximate information on fq” for q = 0. 

The functions were solved mostly by steps of q = 0.2. Since the values 
with q-interval 0.1 must be known for the successive calculations, the 
values lying between were interpolated by means of a Taylor series. For 
the transfer boundary layer the calculations for the cr-value of the air 
(0.7) were performed in the two-dimensional symmetrical case because 
experimental results for the heat transfer of a circular cylinder in air 
exist (see, for instance, the compilation by Rroujillne (ref. l) . For 
(F 0 "\ 0 °he may obtain values from a table given by Goldstein and cal- 
culated by Squire (ref. 2) also in the case of other a- values. Squire 
indicated an analogous expression for the heat transfer at the stagnation 
point. For the rotationally symmetrical case calculations have been 

carried out for a = ^ because the only experimental result for the 

transfer distribution has been found for the evaporation of naphthalene 
spheres (ref. 11), naphthalene has this value of a. If one wants to 
calculate the higher terms of the power-series development in x for a 
special case, one may combine the separate functions in a single term in 
order to save work expenditure; but the generality of the solution is lost 
thereby. This has been done for the boundary layer of the body of revolution 
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for o = — . The parenthesis in the defining equation for Gu, 33j, 

0.395 

etc . , has been combined into a single function and calculated for 

the sphere for the pressure distribution of Page. (See below.) 

The tables of the calculated functions are printed at the end of 
the report with the exception of Idle higher ones for rotationally sym- 
metric boundary layer. Here one has for n = 0 |and or = 

Gg' = 0.3186, E 2* = -0.1005, and %’ = -0.2118. The error of the 

tabulated functions which will be. discussed in more detail later is at 
most a few units in the last digit. 


From the tables one obtains for air, for the pressure distribution 
measured by Hiemenz (ref. 6) at Re ~ 19000 for the circular cylinder 

U = 3.63l4x - 2.l709x 5 - 1 . 5144-x-* 


= 0.9449 - 0.5l00x 2 - 0.5956x 1 *'. . . 

>/5S 

The quantity x is here dimensionless (the length dimension x 
divided by the diameter D) . Not only in Idle range 0° - 55° where the 
series is to apply exactly (see E, 2) but up to the separation point this 
equation is in good agreement with the compilation of experimental dis- 
tribution curves indicated by Kroujiline (ref. l). The derivations will 
be discussed later. 

For the sphere a qualitative agreement with the values obtained for 
evaporation of naphthalene at higher Re (ref —11) is attained if the 
pressure distribution according to Fage (ref. 17) is used which gives 

U = 3x - 3.4966x 3 + 4.7391x 5 - 5.4181x7 for Re = 157200 (ref. 18). 

One then obtains 


= 1.8615 - 2.1477X 2 + 2.4609X 1 *' . . . 

\/Re 

The deviations depend, among other things, on the fact that Fage's 
pressure distribution is possibly not fulfilled for Reynolds numbers as 
small as those used here. Later on a more exact comparison will be made 
with more recent experimental values obtained by the author at still 
higher Reynolds numbers. 
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E. DETERMINATION OP THE ERRORS OF THE 
BROKEN-OFF POWER-SERIES DEVELOPMENTS 


One may use various methods: (l) The following term is calculated 

and for not-too-large x this Indicates the error. (2) If the value 
of the errors is not required with a very high accuracy, the coefficient 
of the x-terms may be assumed to be of the same order of magnitude. (3) 
Use of a continuation method of the profile. (See following section.) 
(4) In the case of direct differentiation, with respect to y, of the 
dc 

value of — taken from the transfer boundary-layer equation one may 
ox 

obtain, for transition to y = 0, the first derivative of Nu in x 
whereby a continuation step may be taken directly with respect to Nu. 
Later on numerical calculations according to seme of these methods will 
be given. 


STEPWISE DEVELOPMENT OF THE BOUNDARY-LAYER PROFILE 


A. TWO-DIMENSIONAL CASE 
a. Flow Boundary Layer 


Prandtl (ref. 8) in d icated for this case a method which is based 
on the fact that one may obtain from the equations an expression for 

^ containing only u with derivatives for a prescribed x. ^ becomes 
with dimensionless quantities without Reynolds numbers 



When two adjacent profiles (at x - Ax and x) were known, for 
instance, by application of the method of Blasius and Hiemenz, it was 
possible to calculate a third for x + Ax. With the u-values at x - Ax 

the 2Ax ^ values for x were vised. One could then continue in the 
ox 

same manner with the profiles at x and x + Ax. In order to guarantee 
the convergence of the expression, one was not to use the original numer- 
ical profile at x but had to replace it by another which satisfied 
certain requirements. In order to calculate those, u was developed 
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into a power series with respect to y: 


a v y v 

u = For the a v one 

v=l v * 


obtained certain conditions by substitution into the flow equation whereby 
only some of them could be chosen arbitrarily. These latter were deter- 
mined by comparison with the given profile. Gortler (ref. 9) perfected 
the method practically and used it in Hiemenz ' pressure distribution over 
the circular cylinder. In the present report corresponding ideas are . 
used for bodies of revolution and for the transfer boundary layer, and 
the necessary expressions are added and discussed. 


b. Transfer Boundary Layer 
From the basic equation one obtains directly 


bo _ l/l d 2 c v dc\ 

bx ulcr by 2 - by) 


This equation may therefore be used directly for step-by-step 
continuation of the vapor and temperature boundary layer. Conditions 
become here simpler insofar as no integration is necessary. However, 
here also the danger exists that the expression becomes uncertain at the 

wall (because of u occurring in the denominator). Moreover, ^ must 

become identically zero at the wall. In order to satisfy the require- 
ments, one resolves here also the quantity c into a power-series devel- 
opment with respect to y(b v function of x only): 


c = 1 


oo 

■E 


V v 


V=1 vi 


By substitution one obtains 

y a v y v y- b v , y v a v , y v ' fl y— b^" 1 _ lr— byy v ~ 2 

v'. v i (v + i) ; ^ — (v - i): a (v - 2)1 

By comparison of terms of the same degree, one arrives at the 
relation between the b v . For the first nine b v there applies (with 
f = -TJU') 
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*1 

a 

a 

^8 


b 9 


*5 


free; ^2 = = 0 

= 2a-jb^’ - a^'b^ free; -jf~ = 3fb^’ - f 'b^; 

= lOa^crt^" + 5a 1 a 1 , (l - 3a)\' + ja-^ClOa 

= 48a 1 crfb 1 " + 2^8^ '(3 - 7a) - l^a^fajb-,/ 
f'a^l^d - 2)[b 1 | 

= 63 f 2 ab 1 " + 7 ff *(2 - 9cr)‘b 1 * + [f , 2 (35<T - 2 ) 


b 6 = 0 

- 1) - a^a^"(5a + l)Jb^ free 
+ [-2a 1 f ,, (3cr +1)- 15fai"cr + 


- ff "(21a + 2) b-L 


The free coefficients are calculated as before by comparison with, 
the given profile, and the c-values developed in power series are sub- 

8c 

stituted into the above equation for — . 

ox 


B. ROTATIOHALLY SYMMETRICAL CASE 
a. Flow Boundary Layer 


The equations read 

u |ii +v |l = U l j» + ^ 
ox dy 

8{url + = o 

8x 8y 

Here r signifies, as before, for a blunt body of revolution the 
distance between the axis of rotation and the base point of the normal 
to the surface. 

By eliminating ^ one obtains a linear equation of the first order 
ox 

in v with solution 



1 dr 

- -r- uy 

r dx 


v 
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By forming the derivative from v with respect to y and using 
the continuity equation one obtains _ 


+ 1 &E. y 
r dx Sy 


For given u-prof ile , one may therefore use an equation for the 
continuation of the boundary-layer profile which differs from the equa- 
tion for the two-dimensional case only with respect to the last term. 

In order to establish the convergence at the wall, here also a power- 
series development in y becomes necessary 


= cL! 
Sx by 





since 


v 



+ 



y v+1 
(v + l)i 


By substitution Into the basic equations one obtains [with f = -UU f 

-3 


a-^ free; a 2 = f; a^ = 0; 

a 4 = a l a l' “ a l 2 S f" ree j a 5 = 2a-^f 1 - 4a^fg; a^ = 2ff ' - 4f^g 


As before, one determines the free coefficients. 


b. Transfer Boundary Layer - 

With the same expression as above for the two-dimensional case one 
obtains for the b v 
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free ; bg = = 0; 

= 28-^' ” a ]_ , \ “ s-^g free; — = Jfb^ - f , b 1 - fb-^g; bg = 0; 


The practical execution in the last three cases will be discussed 
in a later report together with numerical calculations. The methods of 
continuation discussed yield results the accuracy of -which depends 
exclusively on the work expenditure and is therefore not limited by 
postulating approximation functions . The methods may also be used for 
determining the accuracy of the aforementioned power-series developments 
in x in the case of breaking-off after a certain number of terms at a 
certain point. One then starts the continuation method at an x so 
small that the error is- certainly small, and compares the result then 
obtained at a larger x with the one directly calculated from the power- 
series development in x. 


DEPENDENCE OP THE EVAPORATION AND THE HEAT TRANSFER ON a 

A. GENERALITIES 


Pohlhausen (ref. 1 6 ) has shown for the plane that Nu is approxi- 
mately proportional to the quantity \ya. In the approximate calcula- 
tions of Kroujiline (ref. 1) the same was shown for the circular cylinder. 
TJlsamer (ref. 10) demonstrated that the law may be approximately selected 
from various experimental investigations on the heat transfer of a 
circular cylinder. The author of this report has confirmed the law at 
least approximately in the case of evaporation of drops (ref. 11). 

From the equations of the section on power-series developments in 
x one sees that a can probably not be eliminated from them by simple 
transformations. Thus one cannot expect a relation as simple as the 
aforementioned to apply exactly. For the case where a Is very large 
a n d the transfer boundary layer therefore thin compared to the flow 

boundary layer, the author of this report found the ^/~a - law to be exact. 
In this case the curvature of the velocity profile may be neglected in 
the entire transfer boundary layer, and one may replace u by (u 1 )^ 

and v by ^(v")Qy^ in the general houndary-layer equation, with the 
apostrophes indicating derivatives with respect to y. 
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( U > ) r ,y <*£. + ^ V ' . - >0 y2 bs, = 1 ^§2. 

' U 'CV ^ o y ?hr rr ^ O 


Sy a ^y 2 


Hie variable £ = y -^/o" is introduced for c (not for u and v) 

(u . )rf | + <l2o ?2 |.|| 

Boimdary conditions: £ = Oj c = 1. £ = c = 0. 

Thus one has obtained an equation free of a. For this reason, c 

becomes c = f^x,y 34 hence follows that for large a the quantity 

Nu is proportional to the quantity \f/~a. on the entire surface in the 
boundary layer. That the same lav has been found experimentally also 
for a cr that is not large, is based on the fact that the quantity 

Nu W~* does not vary greatly with a and may therefore be found to be 
approximately constant in a small region. 


B. TWO -DIMEIfSIOML SYMMETRICAL CASE 


Into the equations for Fq, F 2, Gfy., % . . . the following 
functions and variables are introduced: 

1 = * FqCtj) = 9 0 (6)} F 2 (ti) = 

Gl*. = Lie 'l l ~ % = 

( 1 )° 

Boimdary conditions: 

|=0; $q = 1; remaining functions = 0; 

£ = oo; all functions = 0. 


j z r -, 4 - ° - 6^(1); etc 
( 1 )° 
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Taking into account that for large cr the equation i(r 
is valid, one obtains 


|(V’)oy 2 


V 


& »» 
w 2 


V 


V 


-5! 2 V 

-3i 2 ®2' + 121*2 
-3I%V + 24|Ti { _ 

-3g 2 e4 f + 24|e^ 


9t 2 ® 0 1 

151 2 * 0 ’ 
151 2 ® 0 ’ + 



For the pressure distribution according to Hiemenz (see p. 20) 

Nu 

one obtains in the case of a circular cylinder -—j=-f== - 1.2592 - 

'5h x y^e 

0.7583x2. . . ; for the case calculated above a = 0.7 one obtains for 

the corresponding qua nt ity 1.0642 - 0.5744-x 2 - 0.6708x^ .... Here 
x signifies the dimensionless abscissa -which is obtained from the length 
dimension through division by the diameter D. The functions 4 >q and 
*2 are given numerically in table 6. 

C. R0TATI0HALLY SYMMETRICAL CASE 


The folio-wing functions are now introduced: 


1 = n ( f i")o J F 0 (n)=* 0 (i)j g 2 (ti) 



The equations become 


V = -3^o’ 

r 2 ” = -3i 2 r 2 ' + 6ir 2 - 6|2* 0 « 

V - "3I 2 9 2 * + 6|9 2 - 6|%o' - | |2 $0 « 
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As in the previous case, the solution of the first equation is 


«b - 1 



For the pressure distribution according to Fage (ref. 17) 
one obtains for the sphere — = 1.^723 - . . )x^ + 

1 


■a = 


0.395 


one obtains 1.3658 - 


(see p. 20) 
. j for 


DISCUSSION OF APPROXIMATE METHODS 


As has been mentioned above, Pohlhausen (ref. 3) gave an approxi- 
mate method for the solution of the boundary-layer equation for the flow 
about a circular cylinder. Tomotika (ref. 18) applied this method to 
the sphere. Kroujouline used a similar method for the transfer for the 
cylinder, applying a broken -off series development in y which was 
determined with utilization of the integral condition of the transfer 
boundary layer. For the flow boundary layer he used a parabolic profile 
whereby the agreement may be assumed to be bad particularly in the case 
of pressure increase. Probably better approximations could have been 
obtained with the us.e of polynominals of the fourth degree. These 
statements are valid only when the transfer boundary layer is thinner 
than the flow boundary layer. Here a brief description is given 
concerning some considerations of the author of this report concerning 
a body of revolution, for various relative magnitudes of the two layers. 

The integral condition formerly not set up for bodies of revolution 
becomes (see p. 4) 


1 d_ 
r dx 



ucdy 



which may be derived, for instance, by integration of the original 
equation. Here 5 is the thickness of the transition boundary layer. 
Using dimensionless quantities without Reynolds numbers only, A Is 

replaced by =k This is assumed below. 
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If for the two profiles the definitions 


S = a l £[ + & 2\i~ 




and 


c = 1 - 2 E + 2 FT - R 

6 \S ) \ 8/ 

are used, one obtains for 8 < 8]_ the equation 


2r = d_ 
8^_za dx 


r 

/ 2 3,4 ' 

___ M a ? z 3 a, z &hZ 

rUSn M= — + -f — + — 3 — + -ft— 

x \15 42 280 18' 


3 

180/ 


Here S^_ is the thickness of the flow boundary layer and z = 


8 

5 l‘ 


For 8 > S]_, the integration is performed, with use of the integral 
condition, first from 0 to 8]_ and then from 8-j_ to 8. Result 


5,zct dx • X| 


2* + P i+fa+fi+fit.J _ 2Q + fa + tl + 5t . 
10 \2 3 4 5 7 z \3 4 5 6 




z3 \ 5 


8 


4 \6 


8 


The two equations have the same form when the parentheses after 
8^ are denoted, for instance, by the letter P. From the first of the 

two equations one sees that, for a cr so large and a z therefore so 
small that only the first term of the parenthesis must be considered. 


this 

zr°- 


z is, for a given x, inversely proportional to the quantity 


Since the Nusselt number Nu equals 




Hu is, for a 


large cr, proportional to fa also according to this approximate 


theory. 


Since r, U, &y , and 8^ are known functions of x, we have in 
any case an equation of the first order with z and x which can be 
solved with customary methods (for instance, with the isocline method 
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or according to Runge and Kutta ) . The only boundary condition required 
for this is the z-value at x = 0. This value is calculated from the 
equation AzaP = 1 where,. A = U'5]_^ is identical with the parameter A 
used by Pohlhausen and Tomotiha. For the sphere where U’ =5 and 
A = 4.716, one obtains in the proximity of the stagnation point 


^ = 0.8759 z^ - 0.2648z^ + 0.0l809z^ + 0.0056lz^ or, respectively, 
1 = 1.4l48z 2 - 1.2295Z + 0.5052 - + ._2i°l§2 


For a given z and therefore also given Nu one may easily 
calculate the corresponding a. For 


Nu 


z - 0.0 0.1 0.4 0-7 1.1 1.6 2.0 5.0 4.0, one obtains 

Re 3/a = 1.526 1.511 1.464 l.4l8 1.356 1.284 1.232 1.128 
1.049 and 4 

l/ 5/a = 0 0.095 O.367 0.622 0.955 1.288 1.545 2.121 2.631 


The quantity Nu j \/Re \fa is therefore, for a large a, almost 
constant and varies in the proximity of the stagnation point about 

linearly with 1 j 


Tie reason for choosing, above, z instead of ,5 as the dependent 
variable was that z probably varies little with x (compare Kroujouline 
(ref. 1)) and can therefore be calculated exactly more easily. 


In the later more detailed report on the investigations, the numer- 
ical results of this formulation as well as of others will be discussed. 
It was shown that the choice of the profile form had a great effect on 
the result. 


SUMMARY 


A preliminary report is given of a theoretical investigation of the 
boundary-layer flow for two-dimensional and rotationally symmetrical 
bodies. The evaporation, the heat transfer, and the velocity are cal- 
culated by power-series developments with respect to the meridian length. 
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The coefficient functions which were calculated numerically in some cases 
have been chosen so that the calculation Is valid for all pressure distri- 
butions and body shapes. The methods for determination of the errors in 
breaking off the series are briefly treated. Methods of continuation are 
discussed. It is shown, for large Prandtl numbers, that the Nusselt 
number is exactly proportional to Hie cube root of the Prandtl number. 
Finally, approximate methods of calculation are discussed. 


EPILOGUE 


I should like to express my deepest gratitude to the director of the 
local institute. Professor Dr. G. Borelius who made working facilities 
available to me and cooperated with me in the friendliest manner during 
the investigation. 

I want to thank also my teacher, the director of the Physical 
Institute of the University Lund, Professor Dr. J. Koch where I carried 
out the previous investigations, for his lasting interest which is so 
very valuable to me. 

To Professor Dr. T. Gustafson, Lund, I am obliged for valuable 
discussions. 


I also owe a debt of gratitude to Kungl. Fysiograf iska Sallskapet 
i Lund which gave me financial support by hiring an assistant for the 
numerical calculations. 

Stockholm., Tekniska Hogskolans Fysiska Ihstitut, February 19^0. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 



32 


MCA OM 1432 


REFERENCES 

1. Kroujiline: Tech. Phys. U. S. S. R. 3 , 183 und 311, 1936. 

2. Squire, not published: Ref. Goldstein, Modern Developments In Fluid 

Dynamics II, 631, 1938, Oxford. 

3. Pohlhausen, K. : Z. S. angew. Math. u. Mech. 1, 252, 1921. 

4. Von Karmen and Millikan: MCA Rep. 504, 1934. 

5. Blasius: Diss. Gottingen, 1907. 

6. Hiemenz: Diss. Gottingen, 1911. 

7. Howarth: Aeron. Res. Comm., Rep &. Mem., Nr. l632,"-1935* 

8. Prandtl: Z. S. angew. Math. u. Mech. 18, 77, 1938. 

9. Gortler: Z. S. angew. Math. u. Mech. 19, 129, 1939* 

10. Ulsamer: Forsch. auf d. Geb. d. Ing. 3, 94, 1932. 

11. Frossling: Gerl. Beitr. zur Geophysik 52, 170, 1938. 

12. Goldstein: s. o. 2., p. 609. 1 

13. Boltze: Diss. Gottingen, 1908. 

14. Geiger-Scheel: Handb. d. Ehys. Ill, p. 631. 

15* Hartree: Cambr. Phil. Soc., Proc. 33, 223 , 1937* 

16. Pohlhausen, E.: Z. S. angew. Math. u. Mech. 1, 115, 1921. 

17. Fage: Aeron. Res. Comm., Rep. &Mem., Nr. 1766, 1937* 

18. Tomotika u. Dual: Proc. Phys.-Math. Soc. Jap. 20, 288, 1938. 



NA.CA rn lk 32 


TABLE 1. TWO-DIMENSIOltAL SYMMETRICAL 
FLOW BOUNDARY LAYER 



*5 

*5* 


® 5 

S5> 

_ If 

85 

*5 


V> 11 

*5 

0.0 

0 

0 

0.7244 

0 

0 

0.6348 

0 

0 

0.1192 

.1 

.0035 

.0675 

.6249 

.0030 



.0005 



.2 

.0132 

.1251 

.5286 

.0114 

.1072 

.4402 

.0017 

.0l4l 

.0249 

.3 

. .0282 

.1734 

-4375 

.0242 



.0032 



.4 

.0476 

.2129 

-3539 

.0405 

.1778 

.2717 

.0045 

.0117 

-.0436 

• 5 

.0705 

.2444 

.2780 

.0595 



.0055 



.6 

.0962 

.2688 

.2112 

.0806 

.2184 

.l4o8 

.0057 

-.0010 

-.0783 

•7 

.1240 

.2869 

.1530 

.1030 



.0052 



.8 

-1534 

.2997 

.1037 

.1264 

.2367 

.0483 

.0039 

-.0176 

-.0833 

• 9 

.1838 

.3080 

.0626 

.1502 



.0017 



1.0 

.2149 

-3125 

.0292 

.1742 

.2399 

-.0106 

-.0012 

-.0330 

-.0680 

1.1 

.2462 

.3140 

.0028 

.1981 



-.0049 



1.2 

.2776 

.3132 

-.0173 

.2218 

.2342 

-.0431 

-.0090 

-.0441 

-.0423 

1.3 

.3088 

.3107 

-.0320 

.2449 



-.0136 



1.4 

•3397 

.3070 

-.0420 

.2676 

.2239 

-.0567 

-.0185 

-.0498 

-.0149 

1.5 

.3702 

.3025 

-.0482 

.2897 



-.0236 



1.6 

.4002 

•2974 

-.0513 

.3112 

.2123 

-.0580 

-.0286 

-.0503 

.0088 

1.7 

.4297 

.2923 

-.0518 

.3322 



-.0336 



1.8 

.4587 

.2871 

-.0506 

.3526 

.2012 

-.0522 

-.0384 

-.0468 

.0256 

1.9 

.4871 

.2822 

-.o48o 

.3724 



-.0430 



2.0 

-5151 

.2775 


.3918 

.1916 

-.0432 

-.0472 

-.0406 

•0351 

2.1 

.5426 

.2733 

-.0402 

.4108 



-.0510 



2.2 

.5698 

.2695 

-.0358 

.4293 

.1859 

-.0335 

-.0546 

-.0331 

.0380 

2.3 

.5966 

.2662 

-.0314 

.4476 



-.0577 



2.4 

.6230 

.2632 

-.0271 

.4655 

.1781 

-.0245 

-.0604 

-.0257 

.0361 

2.5 

.6492 

.2607 

-.0230 

.4832 



-.0628 



2.6 

•6752 

.2586 

-.0194 

.5007 

.1740 

-.0171 

-.0649 

-.0189 

.0312 

2.7 

.7010 

.2568 

-.0160 

.5180 



-.0666 



2.8 

.7266 

-2554 

-.0131 

.5352 

.1712 

-.0114 

-.0681 

-.0133 

.0249 

2.9 

-7520 

.2542 

-.0106 

• 5522 



-.0693 



3.0 

-7774 

.2533 

-.0085 

.5692 

.1694 

-.0072 

-.0703 

-.0089 

.0187 

3.1 

.8027 

.2525 

-.0067 

.5861 



-.0711 



3.2 

.8279 

.2519 

-.0052 

.6030 

.1682 

-.0043 

-.0717 

-.0058 

.0132 

3.3 

.8531 

.2515 

-.0041 

.6198 



-.0722 



3 A 

.8782 

.2511 

-.0032 

.6365 

.1676 

-.0026 

-.0726 

-.0036 

.0089 

3.5 

-9033 

.2508 

-.0024 

.6533 



-.0730 



3.6 

.9284 

.2506 

-.0019 

.6700 

.1672 

-.0015 

-.0732 

-.0021 

.0057 

3.7 

-9534 

.2504 

-.0014 

.6867 



-.0734 



3.8 

.9785 

.2503 

-.0011 

.7034 

.1669 

-.0010 

-.0735 

-.0012 

.0036 

3.9 

1.0035 

.2502 

-.0008 

.7201 


i 

-.0736 



A.o 

1.0285 

.2502 

-.0006 

.7368 

.1668 

-.0004 

-.0737 

-.0006 | 

.0022 

4.1 

1-0535 

.2501 

-.0004 

•7535 



-.0738 



4.2 

I.O785 

.2501 

-.0003 

•7701 

.1667 

-.0001 

-.0738 

-.0003 

.0012 

4.3 


.2500 

-.0002 

.7868 



-.0738 



4.4 



-.0001 

.8035 


-.0001 

-.0738 

-.0001 

.0007 

4.5 

4.6 

4.7 

4.8 



-.0000 

.8201 

.8368 


-.0000 

-.0739 

-.0000 

.0003 









.0001 

4.9 










5-0 









.0000 
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TABLE 2. TWO-DIMENSIONAL SYMMETRICAL 
TRANSITION BOUNDARY LAYER 


n 

£ 

1 

H 

-Ftf’ 

f 2 

f 2 ' 

Gk 

<V 

% 

V 

0.0 

0 

0.4959 

0 

-0.1119 

0 

-0.0977 

0 

0.0318 

.2 

• 0991 

.4953 

-.0224 

-.1113 

-.0195 

-.0970 

.0064 

.0320 

.4 

.1979 

.4917 

-.0443 

-.1077 

-.0386 

-.0932 

.0129 

.0338 

.6 

.2954 

.4825 

-.0651 

-.0988 

-.0565 

-.0846 

*0200 

.0378 

.8 

•3904 

.4660 

-.0835 

-.0841 

-.0721 

-.0704 

.0281 

.0433 

1.0 

.4813 

.4416 

-.0983 

-.0638 

-.0843 

-.0517 

.0374 

.0490 

1.2 

.5666 

.4095 

-.1087 

-.0397 

-.0926 

-.0301 

.0476 

.0529 

1.4 

.6447 

.3708 

-.1141 

-.0138 

-.0964 

-.0078 

.0583 

.0530 

1.6 

• 7146 

.3275 

-.1143 

.0111 

-.0958 

.0131" 

.0685 

.0480 

1.8 

• 7755 

.2818 

-.1099 

.0328 

-.0913 

.0307 

.0771 

.0376 

2.0 

.8273 

.2360 

-.1015 

.0496 

-.0838 

.o44o 

.0832 

.0228 

2.2 

.8701 

.1924 

-.0904 

.0605 

-.0741 

.0521 

.0861 

.0054 

2.4 

.9043 

.1526 

-.0777 

.0656 

-.0633 

.0554 

.0854 

-.0123 

2.6 

.9315 

.1177 

-.0645 

.0654 

-.0522 

.0544 

.0813 

-.0281 

2.8 

.9520 

.0883 

-.0518 

.0610 

-.0417 

.0502 

.0744 

-,o4o4 

3.0 

.9672 

.0644 

-.0403 

.0540 

-0323 

.0440 

.0654 

-.0481 

3.2 

.9781 

.0457 

-.0303 

.0454 

-.0242 

.0368 

.0555 

-.0511 

3.4 

.9858 

.0315 

-.0221 

.0366 

-.0176 

.0295 

.0453 

-.0500 

3.6 

.9910 

.0211 

-.0156 

.0283 

-.0124 

.0227 

.0356 

-.0458 

3-8 

.9944 

.0138 

-.0107 

.0211 

-.0084 

.0168 

.0271 

-.0396 

4.0 

.9966 

.0088 

-.0071 

.0151 

-.0056 

.0119 

.0199 

-.0325 

4.2 

.9980 

.0054 

-.0046 

.0104 

-.0036 

.0082 

• 0l4l 

-.0254 

4.4 

.9989 

.0032 

-.0029 

.0070 

-.0022 

•0055: 

.0097 

-.0191 


.9994 

.0019 

-.0018 

.0045 

-.0014 

• 0035 

.0064 

-.0137 

4.8 

.9996 

.0011 

-.0010 

.0028 

-.0008 

.0021 

.0041 

-.0095 

5-0 

• 9998 

.0006 

-.0006 

.0017 

-.0005 

.0013 

.0026 

| -.0063 

5.2 

.9999 

.0003 

-.0003 

.0010 

-.0003 

.0008 

.0015 

-.0040 

5.4 

1.0000 

.0002 

-.0002 

.0006 

-.0002 

.0005 

.0009 

-.0025 

5.6 


.0001 

-.0001 

.0003 

-.0001 

.0003 

.0005 

-.0015 

5.8 


.0000 

-.0000 

.0002 

-.0000 

.0002 

.0003 

-.0009 

6.0 




.0001 


.0000 

.0001 

-.0005 

6.2 




.0000 



.0001 

-.0002 

6.4 







.0000 

-.0001 

6.6 








-.0001 

6.8 








-.0000 
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TABLE 3- ROTATTONALLT SXMffiOKDCAL PLOW 
BOUNDARY LAYER 


n 

fl 

fl' 

fl" 

s 3 

63’ 

_ 11 

83 

*3 

*3' 

V, 11 

*3 

0.0 

0 

0 

0.9277 

0 

0 

1.0475 

0 

0 

0.0448 

.i 

.0046 

.0903 

.8777 

.0051 

.0998 

.9477 

.0002 

.0044 

.0448 

.2 

.0179 

.1755 

.8277 

.0196 

.1896 

.8488 

.0009 

.0090 

.0444 

.3 

.0395 

.2558 

.7778 

.0427 

.2696 

.7517 

.0020 

.0133 

.0434 

.4 

.0689 

.3311 

.7282 

.0732 

.3400 

.6574 

.0036 

.0176 

.0416 

.5 

.1056 

.4014 

.6788 

.1104 

.4012 

.5666 

.0055 

.0217 

.0391 

.6 

.1490 

.4669 

.6300 

.1532 

.4535 

.4802 

.0079 

.0254 

.0356 

• 7 

.1988 

.5275 

.5819 

.2008 

.4974 

.3986 

.0106 

.0288 

.0314 

.8 

.2544 

.5833 

• 5348 

.2524 

• 5334 

.3227 

.0136 

.0316 

.0265 

■ 9 

.3153 

.6345 

.4888 

.3072 

.5621 

.2528 

.0169 

.0340 

.0210 

1.0 

.3811 

.6811 

.4443 

.3646 

.5842 

.1895 

.0204 

.0358 

.0152 

1.1 

.4514 

.7234 

.4014 

.4239 

.6002 

.1328 

.0241 

.0370 

.0091 

1.2 

.5256 

.7614 

.3604 

.4845 

.6110 

.0832 

.0278 

•0377 

.0032 

1.3 

.6035 

•7954 

.3215 

.5459 

.6171 

.0403 

.0316 

.0377 

-.0026 

lA 

.6846 

.8258 

.2850 

.6078 

.6193 

.0044 

.0353 

.0372 

-.0080 

1.5 

.7686 

.8526 

.2508 

.6696 

.6182 

-.0251 

.0390 

.0361 

-.0127 

1.6 

.8550 

.8761 

.2192 

•7333 

.6144 

-.0483 

.0425 

.0346 

-.0168 

1.7 

• 9437 

.8966 

.1901 

.7925 

.6087 

-.0657 

.0459 

.0327 

-.0202 

1.8 

1.0342 

.9142 

.1637 

•8530 

.6015 

-.0780 

.0491 

.0306 

-.0228 

1.9 

1.1264 

.9294 

.1398 

.9127 

.5932 

-.0857 

.0520 

.0282 

-.0244 

2.0 

1.2200 

.9422 

.1185 

.9716 

.5845 

-.0894 

.0547 

.0258 

-.0254 

2.1 

1.3148 

.9530 

.0996 

1.0296 

-5755 

-.0898 

.0572 

.0233 

-.0256 

2.2 

1.4106 

.9622 

.0831 

I.O867 

.5 666 

-.0876 

.0594 

.0207 

-.0252 

2.3 

1.5072 

.9698 

.0688 

1.1430 

.5580 

-.0834 

.0613 

.0182 

-.0243 

2.4 

1.6045 

.9760 

.0564 

1.1984 

.5500 

-.0776 

.0630 

.0158 

-.0229 

2.5 

1.7024 

.9811 

| .0458 

1.2530 

.5425 

-.0709 

.0645 

.0136 

-.0212 

2.6 

1.8007 

• 9853 

.0370 

1.3069 

.5358 

-.0637 

.0637 

.0116 

-.0193 

2.7 

1.8994 

.9886 

| .0296 

1.3602 

.5298 

-.0563 

.0668 

.0097 

-.0174 

2.8 

1.9984 

• 9912 

.0234 

1.4129 

.5245 

-.0490 

.0677 

.0081 

-.0153 

2.9 

1 2.0977 

• 9932 

.0184 

1.4651 

.5200 

-.0420 

.0684 

.0067 

-.0133 

3.0 

2.1971 

.9949 

.0143 

1.5169 

.5161 

-.0356 

.>0690 

.0054 

-.0114 

3.1 

2.2966 

1 .9962 

.0110 

1.5683 

.5128 

-0297 

.0695 

.0044 

-.0097 

3.2 

2.3963 

.9972 

.0085 

1.6195 

.5102 

-.0245 

. 0699 

.0035 

-.0082 

3.3 

2A961 

• 9979 

.0064 

1.6704 

.5079 

-.0200 

.0702 

.0028 

-.0067 

3 A 

2.5959 

.9985 

.0048 

1.7211 

.5061 

-.0161 

.0705 

.0022 

-.0054 

3.5 

2.6958 

.9989 

.0036 

1.7716 

.5047 

-.0128 

.0706 

.0016 

-.0044 

3.6 

2.7957 

• 9992 

.0026 

1.8220 

.5036 

-.0101 

.0708 

.0013 

-.0035 

3-7 

2.8956 

• 9995 

.0020 

1.8723 

.5027 

-.0078 

.0709 

.0010 

-.0028 

3.8 

2.9956 

• 9996 

.0014 

1.9226 

.5020 

-.0060 

.0710 

.0007 

-.0021 

3-9 

3.0955 

.9997 

.0010 

1.9727 

.5015 

-.0046 

.0711 

.0006 

-.0016 

4.0 

3.1955 

• 9998 

.0007 

2.0229 

.5011 

-.0034 

.0711 

.0004 

-.0013 

4.1 1 

3.2955 

• 9999 

.0005 

2.0730 

.5008 

-.0026 

.0711 

.0002 

-.0009 

4.2 

3.3955 

■ 9999 

.0004 

2.1230 

.5006 

-.0019 

.0712 

.0002 

-.0007 

4.3 

3.4955 1 

■ 9999 

.0003 

2.1731 

.5004 

-.0014 

.0712 

.0002 

-.0006 

4.4 

3.5954 

■ 9999 

.0002 

2.2231 

.5003 

-.0010 

.0712 

.0002 

-.0004 

4.5 

3.6954 

1.0000 

.0001 

2.2731 

.5002 

-.0007 


.0000 

-.0002 

4.6 



.0001 

2.3231 

.5001 

-.0005 



-.0002 

4.7 



.0000 

2.3732 

.5001 

-.0004 



-.0001 

4.8 





.5000 

-.0002 



-.0001 

4.9 






-.0001 



-.0000 

5-0 






-.0001 




5-1 






-.0001 




5-2 






-.0000 
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TABLE 5. ROTATIONALLY SYMMETRICAL 
TRANSFER BOUNDARY LAYER 


a 

-(V)o 

0.5 

•7 

1 

1/0.395 

10 

100 

O.lj-129 

.4705 

.5390 

.7599 

1.2389 

2.7365 


TABLE 6. TWO-DIMENSIONAL SYMMETRICAL TRANSFER 
BOUNDARY LAYER, a LARGE 


1 

1 - $0 


*2 

®2* 

0.0 

0 

1.1198 

0 

-0.4799 

.1 

.1120 

1.1187 

-.0479 

-.4780 

.2 

.2235 

1.1109 

-.0952 

-.4647 

.3 

.3337 

1.0900 

-.i4oi 

-.4293 

A 

.4409 

1.0504- 

-.1801 

-.3637 

• 5 

• 5^30 

.9883 

-.2118 

-.2647 

.6 

.6378 

.9023 

-.2320 

-.1361 

• 7 

• 7228 

•7947 

-.2384 

.0099 

.8 

.7962 

.6711 

-.2301 

.154-1 

.9 

.8567 

.54-02 

-.2083 

.2747 

1.0 

-904-3 

.4-120 

-.1765 

.3530 

1.1 

.9396 

• 2959 

-1395 

•3794 

1.2 

.964-1 

.1989 

-.1023 

.3565 

1.3 

.9801 

.124-5 

-.0693 

.2981 

1.4 

• 9897 

.0720 

-.0432 

.2232 

1.5 

• 9951 

.0383 

-.0246 

.1500 

1.6 

. 9979 

.0186 

-.0128 

.0904 

1.7 

.9992 

.0082 

-.0060 

.0488 

1.8 

.9997 

.0033 

-.0026 

.0236 

1.9 

.9999 

.0012 

-.0010 

.0102 

2.0 

2.1 

2.2 

2.3 

2.4- 

1.0000 

.0004 

.0001 

.0000 

i 

-.0003 

-.0001 

-.0000 

.0040 

.0014 

.0004 

.0001 

.0000 
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